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The phase structure of the 3-d Thirring model 

L. Del Debbio a * (UKQCD Collaboration) 

a Department of Physics, University of Wales Swansea, Swansea SA2 8PP, United Kingdom 

We study the phase structure of the Thirring model in 3-d and find it to be compatible with the existence of 
a non gaussian fixed point of RG. A Finite Size Scaling argument is included in the equation of state in order to 
avoid the assumptions usually needed to extrapolate to the thermodymical limit. 
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1. Introduction 

The Thirring model is a four-fermi theory with 
current-current interaction defined by the follow- 
ing continuum Lagrangian: 



2N 



f 



(1) 



where ip and tp are four-fomponent spinors, m is 
a bare mass, and i runs over Nt fermion species. 

As in the case of the Gross- Neveu model (GN), 
naive power counting suggests that the theory 
is not renormalizablc in the usual perturbative 
expansion in powers of the coupling constant. 
Nonetheless, different approaches might allow a 
continuum limit to be defined for these theories. 
It has been established, using the 1/Nf expan- 
sion, that GN is renormalizable to all orders, 
showing an interesting phase structure with an 
UV fixed point of RG, spontaneous chiral sym- 
metry breaking with dynamical mass generation 
and a spectrum of mesonic bound states. 

For the Thirring model the situation is less 
clear. The 1/Nf expansion turns out to be renor- 
malizable (at least at leading order), but neither 
chiral symmetry breaking nor charge renormaliza- 
tion are observed (the /3-function vanishes identi- 
cally for all values of the bare coupling). On the 
other hand chiral symmetry breaking and dynam- 
ical mass generation are expected at some critical 
value of the coupling g c (Nf), from the solution 
of Schwinger-Dyson equations. However, in or- 
der to solve these equations, one has to use some 
Ansatz for the 1-PI vertex function. According 
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to the chosen Ansatz, the functional dependence 
of g c on Nf turns out to be different. This incon- 
sistent picture can be clarified by a lattice study 
of the model. 

The existence of a non-gaussian fixed point, al- 
lowing a strongly interacting theory with dynam- 
ical mass generation to be defined, would be of 
great theoretical interest. It has a phenomeno- 
logical relevance for model-building beyond the 
Standard Model and it can describe high-T c su- 
perconductivity and the IR behaviour of QED%. 
Besides, due to its relative simplicity, the 3-d 
model is an ideal laboratory to compare differ- 
ent non-perturbative approaches and to test al- 
gorithms for simulating dynamical fcrmions. 

2. Lattice formulation 

We simulate the theory using N species of stag- 
gered fermions and a non-compact action: 



S = 



2 
h.c 



(2) 



For N = 1 this is equivalent to the compact 
formulation where the interaction term becomes 
x(x) e^jp(iA fl (x))x{x + /i) and the mass term for 
the vector field is no longer needed. Both lattice 
formulations have their drawbacks. For N > 1, 
the compact version generates higher order inter- 
actions (six fermions vertices) when the auxiliary 
field is integrated out. The non-compact one has 
an additive charge renormalization, preventing to 
perform numerical simulations at arbitrary strong 
coupling (see p| for details). 
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3. Data analysis: the Equation of State 

The data are obtained, as usual, from numer- 
ical simulations at finite values of the bare mass 
and on finite lattices. In order to study dynamical 
mass generation, we need to consider the m — > 
limit, while we also have to take into account the 
relevance of finite size effects. The only assump- 
tion we are going to make in analyzing our data 
is that the RG flow for the Thirring model is 
similar to the one observed for GN: we therefore 
expect to have a critical coupling g c (N), corre- 
sponding to an UV fixed point, where a chiral 
symmetry breaking phase transition occurs, sep- 
arating a massless phase at weak coupling from 
a massive one at strong coupling. Being below 
the upper critical dimensions (d = 4), one can 
solve the RGE without worrying about logarith- 
mic corrections to scaling. 

An equation of state can then be derived using 
the solution of the RG equations for the exter- 
nal magnetic field at given magnetization H. In 
our case, the external, symmetry-breaking field, 
is the bare mass m, and the magnetization is the 
order parameter monitoring the symmetry break- 
ing, i.e. the chiral condensate. In terms of the 
scaled variables, we obtain: 

m(U)- s = T x (t^V)" 17 ") (3) 

where t = 1 /<? 2 — l/g 2 and T\ is a universal scaling 
function. 

In order to take into account finite size effects, 
we introduce the volume of the lattice as a scal- 
ing field of dimension d. The scaling function de- 
pends now on two independent combinations of 
rescaled variables: 

m(^)- s = T 2 (t(W)~ 1,P ,L~ 1,v {iN>)~ llf> ) (4) 

where L is the linear size of the lattice. 

In the limits L — > oo, t = or m = 0, Eq. (Q) 
reproduces the usual definitions of the critical ex- 
ponents v, 5 and [3. Neglecting the finite-size cor- 
rection and performing a Taylor expansion of T\ 
for small i, one obtains a five parameter equation 
to which the lattice data at fixed lattice volume 
can be fitted, allowing the critical exponents in- 
troduced above to be determined: 

m = B{tpip) + At^) 8 - 1 '? (5) 
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Figure 1. Fisher plot for N = 1, from data at /3 = 
1.6(A), 1.8(<), 2.0(V), 2.2(>), 2.4(+), 2.6(x), on 
a 12 3 lattice. 



Our results for JV = 1, 2 are reported in Tab (0), 
where Fit I corresponds to Eq. (||), while for Fit 
II we set 5 — X//3 = 1, as suggested in ||]. It 
is intersting to remark that, within errors, the 
hypothesis underlying Fit II is compatible with 
the result obtained with (3 as a free parameter. 
Since Eq. (|^) is obtained from a Taylor expansion, 
we expect it to be valid only for a limited range 
of t values. The number of points included in 
our fits has been chosen in order to minimize the 
X 2 /dof. The good quality of the fit can be seen 
as an evidence in favour of the description of the 
phase structure of the Thirring model in terms of 
a non-gaussian, UV stable, fixed point. Moreover, 
the critical exponents are significantly different 
from their mean-field values; i.e. the continuum 
theory defined by the lattice model is not trivial. 

If the critical behaviour is described by mean- 
field theory, the equation of state becomes: 

— 771 

(^) 2 = Kl7T - + K 2 t (6) 

so that (■ipip) 2 is a linear function of the ratio 
mj (tptp) and a positive value of the intercept cor- 
responds to a non-vanshing value of the chiral 
condensate for m — 0, while the intercept will be 
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Table 1 

Results from fit 



Parameter 


Fit 1 


Fit 11 


Fit III 


Mean Field 


N = 1 1/gi 


2.03(9) 


1.94(4) 


2.03(3) 




S 


2.32(23) 


2.68(16) 


2.65(7) 


3 




0.71(9) 






1/2 


V 






0.80(8) 


1/2 


A 


0.32(5) 


0.37(1) 


0.33(9) 




B 


1.92(43) 


2.86(53) 


2.33(19) 




X 2 /dof 


2.4 


2.1 


1.5 




N = 2 l/.g c 2 


0.63(1) 


0.66(1) 






5 


3.67(28) 


3.43(19) 




3 




0.38(4) 






1/2 


A 


0.78(5) 


0.73(2) 






B 


7.9(2.8) 


6.4(1.5) 






X 2 /dof 


3.1 


2.0 







exactly zero at the critical coupling. This plot, 
known as Fisher plot, is shown in Fig. ([!]) for 
N = 1, where a clear evidence for chiral symme- 
try breaking can be seen. 

Finally, one can include finite size effects into 
the analysis, by expanding JF 2 in Eq. (||). After 
setting 6 — 1/(3 = 1 in order to reduce the num- 
ber of parameters, an equation to fit data from 
different lattice sizes can be written as: 

m = A(^) s + B(t + CL- 1 / U )(ijip) (7) 

The results of the fit for N = 1 are reported in 
Tab. (0) as Fit III. They include data for differ- 
ent values of t and m, from 8 3 , 12 3 and 16 3 lat- 
tices. The critical exponent r\ is defined in terms 
of (5 and the space-time dimension d using hy- 
perscaling, while (3 has been fixed by the above 
constraint. Using the results from Fit III, /?, v 
and r\ verify, within errors, the scaling relation: 

P=\v{d-2 + n) (8) 

confirming the fact that our data can be described 
by a UV fixed point in the RG flow of the model. 

4. Conclusions 

Our results can be summarized as follows. We 
have found a numerical evidence for a chiral sym- 
metry breaking phase transition in the Thirring 



model with TV = 1,2, showing a breakdown of 
the leading order 1/N prediction for small N. A 
more careful analysis of the N = 3 case could also 
shed some light on the different scenarios pre- 
dicted by the Schwinger-Dyson approach. Our 
data are consistent with the existence of a non- 
gaussian fixed point, corresponding to this phase 
transition. By combining the equation of state 
with finite size scaling analysis, we are able to 
determine the critical exponents. These indicate 
that the continuum theories defined at the critical 
points are not mean-field. In order to clarify this 
structure, we are now studying the spectrum of 
the theory in the different phases, the renormal- 
ized coupling and the lines of constant physics. 
A detailed description of our results will appear 
soon 0. 
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